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Abstract 

We use the Dirac equation to study the mass spectrum of mesons 
with heavy-hght quark combinations. First we study the Dirac equation 
with sphericaUy symmetry and funnel potential, and apply them on the 
hydrogen-like atom problem to check the correctness of our numerical 
program. Then we test the parameters in Olsson's paper [T]. We show 
that Olsson's parameters are good in fitting the averaged central mass, but 
fail to get correct energy fine splitting. Finally we fit the mass spectrum 
data of D, Ds, B and Bs mesons with our parameters by solve the Dirac 
equation and funnel potential, calculate the energy splitting of the S and 
P states. Our parameters can fit the mass and fine splitting with errors 
in less than IMeV . 

PACS numbers: 12.40.Yx, 14.40.-n, 12.39.-x, 03.65.Pm 
Keywords; Meson, Quark, Mass Spectrum, Dirac Equation. 

1 Meson Mass Spectrum Question 

In the standard model, a meson is composed of a quark and an anti-quark, bound 
together by the strong interaction. Through the studying of the mass spectrum 
of mesons, we can demonstrate the correctness of the quantum fields theory, 
and predict the particle's mass that has not been found yet in the experiments. 

Many articles had studied the meson mass spectrum problems and got many 
very good results [I] [1] [3] [1]. Most authors used the Schrodinger equation to 
solve the problem. Their treatment is accurate enough for the heavy mesons, 
which are composed by two heavy quarks and moving slowly, thus may be 
treated non-relativistically. Let's take an estimate. If we think the mass of a 
meson is the combination of the total mass, kinetic and potential energies of the 
two composition quarks, the binding energy (kinetic -|- potential) is calculated 
and listed in Table [T] The constituent mass values we used are: rriu = 0.30 



•email: pan_b@sina.com 

^ Member of the BaBar Collaboration at slac.stanford.edu since 2004. 



1 



GeV, TOd = 0.30 GeV, = 0.45 GeV, = 1.3 GeV, rrib = 4.5 GeV and 
m-t = 180 GeV. The "ratio" column is defined as 

binding energy 

Ratio = — — (1) 

light quark s mass 

Compared to the light quark's mass in the J/V', T and mesons, their kinetic 
and potential energies are not very large. But in the Z?, Dg, B and Bg mesons, 
the ratio is very large. So the light quarks in these mesons, £), Ds, B and 
Bg^ are moving with relativistic energies. Thus requires us to use relativistic 
equation, the Dirac equation, to solve the spectrum problem. 



Table 1: The composition quarks and mass of some mesons fT] [21 [SI. 



Meson's 


Meson 
mass 
(GeV) 


Composition 


Binding 
energy 
(GeV) 


Ratio 


name 


quarks 


D 


1.975 


cu 
cd 


0.37 


1.23 


Ds 


2.075 


cs 


0.32 


0.75 


B 


5.314 


ub 
db 


0.37 


1.23 


Bs 


5.410 


sb 


0.32 


0.75 




3.097 


cc 


0.40 


0.28 


T 


9.464 


bb 


0.40 


0.09 


Be 


6.264 


cb 


0.40 


0.28 



We should point out that the quark's mass we used here is the constituent 
quark mass, which is the quark's current mass plus the mass of the gluon fields 
and sea-quarks, as the effective quark mass of the valence quark. The current 
quark mass means the mass of a quark itself only. Values of the composition 
mass and current mass differ greatly. For example, proton's mass is about 0.938 
GeV. The rest current masses of its three valence quarks are only about 0.011 
GeV each. But we can treat the mass of each up or down quark with constituent 
mass as large as 0.30 GeV. 

In a charm meson, there are two quarks, with rric ^ rriu^d- So we may 
simplify the problem by treat the heavier quark's mass as infinite. Then the 
problem is reduced to the light quark moving in the funnel potential that created 
by the heavier quark. Solving the Dirac equation with the funnel potential, we 
can get the eigenenergy of the light quark in the funnel potential. Because of the 
Dirac equation includes the light quark's spin-orbit interaction, the eigenenergy 
will refiect the splitting of the S and P states. The influence of the heavier 
quark is treated by adding its mass in the spin dependant force. The total mass 
of the meson can be expressed as: 

Mg^g^ = + nig^ + E + A, (2) 

that's the sum of the masses of the two quarks rug-^ and rriq^ , eigenenergy of the 
system E, and the hyper-fine energy splitting A. 
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M. G. Olsson [T] shows that the Dirac equation can be used to get a good 
fitting in the average mass spectrum. But he did not calculate the fine structure 
splitting. 

D. Ebert, V. O. Galkin and R. N. Faustov |2] use Schrodinger equation with 
relativistic potential to solve the mass spectrum and calculate the fine splitting. 
They point out that in the spin dependent potentials, we should use eigenlnergy ' 
instead of use in Eichten and Quigg's paper [3^ to get energy splitting. 

In this paper, we will solve the Dirac equation numerically to fit the meson's 
mass spectrum of the heavy- light quarks combination system, such as D, Ds, 
B and i?s, and calculate the spin-orbit energy hyper fine splitting. Our result 
will be compared to the experimental data. 



2 Dirac Equation with Spherically Symmetry 



The Dirac equation for free particle is 

d 



(7m^ + w)V' = 0. 



(3) 



If the potential has spherical symmetry, and can be written as a combination 
of Lorentz scalar part Vs{r) and vector Coulomb potential Vy{r), then the Dirac 
equation can be written as [B] [7] [5]: 



d 



ia, — + /3{m + Vs{r)) + Vy{r) 



with the solution has the form like 



Fir) 
G{r) 



r{r)y%Af) 



(4) 



(5) 



in which yjj^ is a two component generalization of the spherical harmonic func- 
tions Yij^, 



k 

The number k is 



2k + I 



k = ±[j 



2k + 



in which the ± sign is 



if l^j 



if 1^3-^ fc = -(j + i) + 

Thus the coupled equations for the radial functions are 

dg{r) 1 — /c 



E-m-Vs{r)^V,ir) f{r) = — 



dr 



(6) 
(7) 

(8) 

(9) 
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E + m + Vs{r) - 

Notes that the angular dependent variables have been removed, and only 
two radial coordinate dependent functions / and g are left in the equation. By 
solving the equations © and PU|) . we will get the energies and eigenfunctions 
of a relativistic particle that moves in the central potential. 



dr 



-fir)- 



(10) 



3 Potential inside Meson 

3.1 QCD One Gluon Exchange Coulomb-like Potential 

Quantum Chromodynamics (QCD) describes the strong interaction among quarks 
and gluons. The Lagrangian density in QCD can be written as: 

^ = - J^;.^r + E 97(i p - mf)qf, (11) 



/=i 



in which 



= d^.-igsTaA% 

where qf is the quark field, gs is the strong interaction coupling strength, Ta 
are the generators of the color SU{3) group, Sabc is the structure constant of 
the su(3) Lie algebra, and A^ is the SU{3) color gauge field. 

Effective potential is used to study the bound states of the quarks in mesons. 
The one gluon exchange Coulomb-like potential jlOj between two quarks is 

VAr) = -^^, (12) 
6 r 

in which 

a. = |i. (13) 



3.2 Funnel Potential 

At large distances, there should exist a confine potential that describes the color 
confinement. Unfortunately, till now, the color confinement can not be derived 
from the QCD first principle. So we may add it in "by hand". The commonly 
used Cornell confinement potential [4] [11] [12] is 

Vs{r) = Vconfmeir) = QT, (14) 

which will make the combination of the scalar Vs (r) and vector Vy (r) potential 
has a funnel shape. The spectrum obtained with the funnel potential is in good 
agreement with experimental data for the light and heavy mesons [1] [2] [2]- 
Typical values for the parameters are a « 0.2 GeV"^, and as ~ 0.2 ^ 0.3. 



4 



4 Spin Dependant Hyper-fine Energy Splitting 



In QED, the relativistic Dirac equation describes two or more massive spin 1/2 
particles interacting electromagnetically. The perturbation QED formula of the 
fine splitting includes spin-orbit interaction (L-S), spin-spin interaction (si-Sj), 
and Breit spin-spin tensor interaction (si ■ Sj — LlLliil^Ej.ljjl'j terms. 

ij 

Within the framework of QCD, there are similar potentials. If the Schrodinger 
equation is used to solve the mass spectrum problem, the spin-dependent hyper- 
fine energy splitting to the first order of as can be written as [3j: 

4 

A = Vr, (15) 



k=l 



with 



2mf 2mj 



and 



rp {L- S^)~ {L-Sj)~ 

miirij mirrij 

{Si ■ Sj) ~ 

Ts = T3{mi,mj), 

rmmj 

T4 = ^^f4{m,,m,), (16) 
rmmj 



,ldV^ 

Ti[mi,mj) = -{-—} +2T2[mi,mj), 

Um,,m,) = (17) 



in which 



Sij = 4[3(si • n){sj ■ n) - Si ■ Sj]. (18) 

Now we use the Dirac equation to solve the mass spectrum. The energy 
splitting terms in equ. (|16p and (fT7|) should be modified a little bit. In the non- 
relativistic limit, p — mv << nic. If we let Vg = Q and V — Vy, the Dirac 
equation dU and wave function ^ will be reduced to: 



p2 



dv dF n idv - , , , 9. . . 

4m^c^ ar ar Im^c r ar 



(19) 
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Compare to the Schrodinger equation, the non-relativistic hmited Dirac 
equation ([19]) has aheady included the 2mir ' ^ spin-orbit interaction term. 
So when we use the Dirac equation to study energy sphtting, the — (f ^) term 
should be removed from Ti(mi,mj) in equ. (|17p . Since we treat the heavier 
quark's mass as infinite in the Dirac equation, the fracture will be replaced 
by a meaningful limited number: over energy eigenvalue, Now the terms in 
the energy splitting eau. pS)) will be 

EiTij Errij 
Enij ' 

r, ^ IMf,, (20) 



and 



T2[m,,mj) = — (r ), 

T4(m„m,) = (21) 

There is a trick in calculating the ip{0)'^ term in equ. (PT|) . It can be replaced 

by the production of the wave function and a S function: (5(r)^(r)^. The S 
function can be defined as 

<5(r)^ lim -i^e-f^. (22) 



So 



9 
9 







6^0 963^1 Jo 

which can be easily integrated numerically. 

By adding the spin dependant part (|T5t with pO|) and ((2T|) . into the solution 
of equ.® and equ. lfTO)) . we will get the meson spectrum with fine structure 
energy splitting. 
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5 Numerical Results for the Mesons 



We will use the double shooting method and Runge-Kutta 4th method to solve 
the Dirac equation. Because the Dirac equation for the hydrogen-like atom 
has exact analytical solution, we first run our program on the hydrogen-like 
atom case to test whether our code works or not. The test results are shown in 
Appendix A, that we can get up to f accuracy. Next we run with Olsson's 
parameters and the funnel potential to compare to Olsson's [T] results, which 
are shown in Appendix B. 

Finally with the funnel potential 

Vs{r) - ar, (24) 

V,{r) = (25) 
r 

we solve the Dirac equation numerically with our own parameters. Because 
meson has a confining part of potential, the wave function tends to contract to 
the center. So the system will be smaller in scale than the hydrogen-like atom, 
which does not have a confining part. According to our test, we choose the 
boundary condition as 

?/'('') lr=0 (Gey- 1) = 0, (26) 

and 

V'(f')|r>20(Gey-i) = 0. (27) 

The radial part of the wave functions will be normalized with the Simpson's 
integration rule to one, 

poo 

if+g^ydr^l. (28) 



The fitting method is done by inputting two mesons masses as initial values 
to determine the parameters in the expression ([T5l) , equ ([9|) and ([TOl) . Then use 
the trial parameters, input configuration parameters of an unknown meson, to 
get the mass of that unknown meson. 

After many trials, we find out the following set of parameters can fit the 
meson's average mass and splitting very well. 

mud = 0.280 GeV, 

nis = 0.429 GeV, 

rUc = 1.095 GeV, 

rub = 4.435 GeV, 

a = 0.368 GeT/^ 

K 0.400. (29) 

Our fitting results are listed in Tabled in which the particle's experimental 
mass are from the PDG (5j book. Spin averaged mass is calculated by taking | 
(|) of the triplet mass, and j (|) of the singlet mass for the s(p) states [1]. The 
column " Numerical center mass" are the numerical result of the central mass of 
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Table 2: Using Dirac equation to fit the meson spectrum. Particle's experimen- 
tal mass are from the PDG 5 book. Mass values are in M eV. 



States 


2S+1 T 


Spin 

averaged 
mass 


7^ 
K 


Numerical 
center 
mass 


Numerical 
splitting 
mass 


Parameter 
b 


15(1867) 
L) (2010) 


'So 


1S(1974) 


-1 


1975 


1867 
2010 


b=1.07 


L>i(2423) 
15^(2457) 


'Pi 

'P2 


2P(2444) 


-2 


2407 


2423 
2457 


N/A 


£'3(1969) 
i:'*(2110) 


'So 

's. 


1S(2075) 


-1 


2074 


1968 
2109 


b=1.08 


£',1(2535) 
^.j(2573) 


'Pi 

'P2 


2P(2559) 


-2 


2515 


2528 
2570 


N/A 


5(5279) 
B*(5325) 


'So 
'Si 


1S(5314) 


-1 


5314 


5279 
5325 


b=0.87 


5,(5375) 
B;(5421) 


'So 

's. 


1S(5410) 


-1 


5412 


5376 
5422 


b=0.88 



the S and P states. Then we use the fine structure formula (fTS)) to calculate the 
energy fine splitting that are listed in the column "Numerical splitting mass". 
For the P states, by intentionally choosing parameters that let the spin average 
mass does not sit between the ^Pi and 'P2 states, but let the average mass 
lower than both of the ^Pi and 'P2 states, we can get good fittings for their 
splitting. The errors for S states are about 1 MeV, while the P states errors 
are less than 7 MeV. 

We also calculate the average values of r, r^, and r~^, which are listed 
in Table [31 In our model, the wave functions are related to the light quark's 
mass, but not to the heavy quark's mass. 

6 Discussion 

By using our set of parameters (P^)) . we can fit the D, Dg, B and Bg mesons 
spectrum with the errors within ^ 7 MeV . The parameters for the mass of the 
u, d, s, c and b quarks, k are all in reasonable range. We use a = 0.368 GeV^. 
But in other people's paper [2] [3], a « 0.18 ~ 0.20 GeV^, in which they used 
with the Schrodinger equation. The reason may be explained as the differ- 
ence between the Dirac equation and the Schrodinger equation. The one-gluon 
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Table 3: Average value oi r, r^, r ^ and r ^, with values are in GeV^' 



Light quark '^^^^Lj <r > <r^> <r^> <r^> 



u/d 
u/d 



'Pi 

'P2 



1.520 2.811 



2.205 5.423 



0.887 



0.521 



1.412 



0.326 



^Sq 



1.437 



2.521 



0.944 



1.617 



'Pi 

'P2 



2.120 



5.026 



0.543 



0.355 



exchange process plus the confining potential is 



V{p,q) 



-u{p) 



4q!s 47r 
'~3~fc2 



7° + Vconf 



in which 



and 



u{p) 



Ep + m 
2m 



1 

Ep+m 



+ 

u u = \ — . 
V m 



(30) 



(31) 



(32) 



In the non-relativistic limit, assume the exchanging gluon's energy is small, 
then Ep K, Eq. By using E"^ = + rn^, we can get 



4as 47r m 



(33) 



That means when we use confining potential Vconf in the Dirac equation, it is 
equivalent to the confining potential ^Vconf in the Schrodinger equation, 



TO TO 
Vconf = ar > Vconf = "^or. 



(34) 



So the relation between the parameters of "a" in the Dirac equation and the 
Schrodinger equation is 



TO 



Dirac: a <^=^ Schrodinger: — a. 

^ E 



(35) 



Let's take an estimate. For the D meson, using our parameters ([29]), Maverage 
1.975 GeV, to„ — 0.28 GeV, TOc = 1.095 GeV, so the eigenenergy is 



E = Ma 
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Wave functions of the D meson with Our parameters 



CO 

c 
o 



i 

"5 
r 

n 

Q. 
« 

(0 
OC 



1.5 



0.5 



-0.5 





n = l,k = -l 






j = 0.5, I = 






kapa = 0.400 






/ 

g 








. , , . 1 . . . , 1 . . , , 1 . . . . 1 , . . , 1 . . . . 1 , , . . 1 . . . , 1 . , , . 1 . . . 



8 9 10 
r (GeV^) 



Figure 1: Wave functions of the IS* state with our parameters. 



= 1.975 - 1.095- 0.28 
= 0.60 GeV. 

That means the parameters of " a" in the Schrodinger equation is 
™ X 0.368 = 0.184 (GeF)2, 



E 0.60 

which is in the range that people used with the Schrodinger equation. 



(36) 



(37) 
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Wave functions of the D meson with Our parameters 




Figure 2: Wave functions of the 2P state with our parameters. 



Appendix A: Numerical Results for the Hydrogen- 
like Atoms 

The hydrogen-hke atom is defined as a particle moves in the central Coulomb 
potential. We can use its analytical solution results to test the correctness of 
our numerical program. 

For the hydrogen-like atom with a Coulomb central potential 



V 



K 

r 



the exact Dirac solution [7] [5] is 



En,k = m 



(iV+|fc|)2 + 2iV(Vfc2-K2_ |fc|) 



1 - 



where 



n = N+\k\ > 1, -n<k<n. 
We use the following parameters in our numerical code. 

m = 0.3 GeV, 



(38) 



n2 + 2{n - ij + i)) J{j + 1)2 - «2 _ (j- + 1) 



,(39) 



(40) 
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Za = K = 0.579, 

V = (41) 

r 

Because the hydrogen-like atom does not have a confining part of potential, 
the wave functions tend to extend far away from the center. So the system will 
be large in scale. In our numerical program, we should set a large region to 
solve the problem. According to our test, we choose the boundary condition as 

V'(Olr=0(Gey-i) = 0, (42) 

and 

V'('')lr>200 (GeV-i) 0. (43) 

Both of the numerical and analytical results are listed in Table ID It shows 
that our numerical program works very well, that there are no differences among 
the numerical and analytical results of the eigenenergy to the precision of 10~^. 

We also calculate the average value of r, r^, and r~^, which are listed 
in Table [Slin] and [T] Our numerical results agree with the exact Dirac solution's 
expect values. The Schrodinger exact solution's expect values are also listed for 
comparison. In the non-relativistic limit, Dirac equation's / part wave function 
will reduce to the Schrodinger wave function. The difference in our results 
between the Dirac and Schrodinger equations indicate that for a system with 
quark's mass and interaction strength, we should use relativistic Dirac equation, 
instead of the non-relativistic Schrodinger equation. 



Appendix B: Numerical Results of with Olsson's 
Parameters 

The parameters in Olsson's paper [1] are 

mud = 0.300 GeV 
m, = 0.463 Gey 
mc = 1.301 GeT/ 
mb ^ 4.639 Gey 
0.308 GeT/^ 

0.579, (44) 



with the potential 



a 

K 



Vs{r) = ar, (45) 



K 



V^{r) = . (46) 



We solve the Dirac equation with Olsson's parameters p4|) . Our fitting re- 
sults of the meson's mass spectrum for both center average mass and energy 
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Table 4: Exact Dirac solutions and our numerical results of hydrogen-like atom. 



n 


k 


j 




name 


Analytical energy (Gel^) 


Numerical 
result 
) 


1 


-1 


1 

2 





l'S'l/2 


mVl - Za^ =0.24460 


0.24460 


2 


-2 


3 
2 


1 


2-P3/2 


m^l - i(Za)2 = 0.28715 


0.28715 




1 

-1 


1 

2 

1 
2 


1 




2^^1/2 ] 

251/2 J 


\mjl 0.28581 


0.28581 


3 


-3 


5 
2 


2 


3D5/2 


ra^\ - \{ZaY = 0.29436 


0.29436 




2 

-2 


3 
2 

3 
2 


2 
1 


3D3/2 ] 
3P3/2 J 


I m, /l ^^S2^ - 0.29425 

V 5 + 2./4-(Za)2 


0.29425 




1 

-1 


1 

2 

1 
2 


1 



3Pl/2 ] 

35i/2 J 


I m, /l 0.29385 

I V 5+4Vl-(Za)2 


0.29385 



splitting are listed in Table [5] with Olsson's parameters. The results show that 
we can reproduce the spin average center mass values, which are listed in Osson's 
paper. Because we have the parameter 6, which is the width of the (5 function, 
we can get good results for the fine splitting of the S states by adjusting the 
values of h. On the other hand, the calculated values of the splitting for the P 
states are unfortunately not so good, with the errors are around 90 MeV. 
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Table 5: Average value of r and of the hydrogen-like atom. Exact Dirac 
(column "Exact(D)") and Schrodinger (column "Exact(S)") solutions are also 
listed for comparison to our numerical results. 



n 


k 


j 


1 


Numeric 


<r> {GeV- 
Exact(D) 


Exact(S) 


Numeric Exact(D) Exact(S) 


1 


-1 


1 

2 





7.572 


7.572 


8.635 


79.14 


79.14 


99.43 


2 


-2 


3 
2 


1 


27.80 


27.80 


28.79 


932.9 


932.8 


994.3 


2 


1 


1 
2 


1 


24.24 


24.25 


28.79 


729.7 


730.4 


994.3 


2 


-1 


1 
2 





29.99 


30.01 


34.54 


1074. 


1074. 


1392. 


3 


-3 


6 
2 


2 


59.44 


59.47 


60.45 


4043. 


4050. 


4176. 


3 


2 


3 
2 


2 


57.71 


57.73 


60.45 


3832. 


3836. 


4176. 


3 


-2 


3 
2 


1 


69.16 


69.24 


71.96 


5540. 


5558. 


5965. 


3 


1 


1 
2 


1 


64.14 


64.18 


71.96 


4798. 


4806. 


5966. 


3 


-1 


1 
2 





69.87 


69.94 


77.72 


5597. 


5611. 


6861. 
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Table 6: Average values of in GeV of the hydrogen-like atom. Exact Dirac 
(column "Exact(D)") and Schrodinger (column "Exact(S)") solutions are also 

listed for comparison to our numerical results. 

n k j 1 Numeric Exact (D) Exact (S) 

1 -1 i 2.13 X 10"^ 2.13 X 10-1 1.74 x 10"^ 

2 -2 I 1 4.54 X 10-2 4.54 x 10"^ 4.34 x 10"^ 
2 1 i 1 5.59 x 10-2 5.59 x lO'^ 4.34 x 10"^ 

2 -1 i 5.59 X 10-2 5.59 x lO'^ 4.34 x lO'^ 

3 -3 I 2 1.97x10-2 1.97x 10-2 1.93 x 10-2 
3 2 I 2 2.05 X 10-2 2.05 x 10-2 1.93 x 10-2 
3 -2 I 1 2.05 X 10-2 2.05 x 10-2 1.93 x 10-2 
3 1^ 1 2.36 X 10-2 2.36 x 10-2 1.93 x 10-2 
3 -1 i 2.36 x 10-2 2.36 x 10-2 1.93 x 10-2 
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Table 7: Average value of in GeV^ of the hydrogen-like atom. Exact Dirac 

(column "Exact(D)") and Schrodinger (column "Exact(S)") solutions are also 
listed for comparison to our numerical results. 

n k j 1 Numeric Exact (D) Exact (S) 

1 -1 i 1.15x10"^ 1.17x10-1 6.03x10-2 

2 -2 I 1 2.79 X 10-3 2.79 x 10-^ 2.51 x 10-^ 
2 1 i 1 5.82 x 10-3 5.84 x 10-^ 2.51 x 10-^ 

2 -1 i 1.83 X 10-2 1.87 X 10-^ 7.54 x 10-^ 

3 -3 I 2 1.97x10-2 1.97x10-2 1.93x10-2 
3 2 I 2 5.14 X 10-4 5.14 x 10-"^ 4.47 x 10-" 
3 -2 I 1 8.67 X 10-^ 8.66 x 10-* 7.45 x 10-* 
3 1^1 1.79 X 10-3 1.80 X 10-3 7.45 x 10-* 
3 -1 i 5.44 X 10-3 5.56 x 10^3 2.23 x 10-3 
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Dirac wave functions of the hydrogen like atom 
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Figure 3: Numerical results of the Dirac wave function of the hydrogen-like 
atom's 15*1/2 state. 
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Figure 4: Numerical results of the Dirac wave function of the hydrogen-like 
atom's 2Pi/2 state. 
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Dirac wave functions of the hydrogen like atom 
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Figure 5: Numerical results of the Dirac wave function of the hydrogen-like 
atom's 25*1/2 state. 
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Figure 6: Numerical results of the Dirac wave function of the hydrogen-like 
atom's 2P3/2 state. 
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Dirac wave functions of the hydrogen like atom 
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Figure 7: Numerical results of the Dirac wave function of the hydrogen-like 
atom's 351/2 state. 
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Figure 8: Numerical results of the Dirac wave function of the hydrogen-like 
atom's 3Pi/2 state. 
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Dirac wave functions of the hydrogen like atom 




Figure 9: Numerical results of the Dirac wave function of the hydrogen-like 
atom's 3P3/2 state. 



Dirac wave functions of the hydrogen Wke atom 



0.0035 



n = 3.k = 2 
j = 1.5, I = 2 
m = 0.3 GeV 
kapa = 0.579 
-f 




100 120 140 160 180 200 

r (GeV^) 



Figure 10: Numerical results of the Dirac wave function of the hydrogen- like 
atom's 3D3/2 state. 
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Dirac wave functions of the hydrogen like atom 




Figure 11: Numerical results of the Dirac wave function of the hydrogen- like 
atom's 3I?5/2 state. 



Table 8: Using Dirac equation to fit the meson spectrum with Olsson's param- 
eters (j44l) [1 . The values of mass are in MeV . 


States 




Spin 

averaged 
mass 


k 


Numerical 
center 
mass 


Numerical 
splitting 
mass 


Parameter 
b 


15(1867) 
i:>*(2010) 




1S'(1974) 


-1 


1975 


1867 
2010 


b==1.64 


£'1(2423) 
15^(2457) 


'Pi 

'P2 


2F(2444) 


-2 


2444 


2535 
2428 


N/A 




'So 


1S'(2075) 


-1 


2074 


1968 
2109 


b=1.61 


Dsi{2535) 


'Pi 

'P2 


2P(2559) 


-2 


2559 


2537 
2660 


N/A 


5(5279) 
B*(5325) 


'So 


1S'(5314) 


-1 


5314 


5279 
5325 


b=1.43 


5^(5375) 
B,*(5421) 


'So 


1S'(5410) 


-1 


5412 


5376 
5422 


b=1.43 
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Table 9: Average value oi r, r^, r ^ and r ^, with values are in GeV" , by using 
Olsson's parameters. 

Light quark ^^'+^Lj <r > < > < r^^ > < r"^ > 
To"^ 

u/d 3";" 1.482 2.736 0.956 1.899 
1 p 

u/d 3„ 2.286 5.867 0.507 0.313 

^2 

"'"'So 



1.369 2.353 1.046 2.306 



2.177 5.333 0.534 0.349 




Figure 12: Wave functions of the 15 state with Olsson's parameters. 
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Wave functions of the D meson with Olsson's parameters 




Figure 13: Wave functions of the 2P state with Olsson's parameters. 
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